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Abstract 

We compute the energy of a static hybrid, i.e. of a hybrid quarkonium with static quark and 
antiquark, at short distances in Z) = 4, 3 dimensions. The soft contribution to this energy is the 
static potential of a color octet quark-antiquark pair at short distances, which is known at two 
loops for arbitrary D. We have checked this expression employing thermal field theory methods. 
Using the effective field theory pNRQCD we calculate the ultrasoft contributions to the hybrid 
(and singlet) static energy at the two-loop level. We then present new results for the static hybrid 
energy/potential and the hybrid decay width in three and four dimensions. Finally we comment 
on the meaning of the perturbative results in two space-time dimensions, where the hybrid does 
not exist. 

PACS numbers: 12.38. Cy, 12.38.Bx, ll.lO.Kk, ll.lO.Hi 
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I. INTRODUCTION 

In a recent paper [l| we have studied the potential and energy of a static color-singlet 
quark- ant iquark state at short distances in D space-time dimensions (though with special 
emphasis on the three-dimensional case). It is the aim of this paper to perform a similar 
analysis for the static hybrid energy and the associated octet potential. 

The energies of static hybrid systems are very interesting quantities. They can teach 
us a lot about the dynamics of QCD, and can be potentially relevant for the theoretical 
description of physical hybrids made of heavy quarks. Their behavior at long distances 
might contain information on the detailed dynamics responsible for confinement, whether 
it is due to strings and, if so, of which kind. Most relevant for us, however, is that the 
hybrid energy at short distances approximately equals the sum of the static octet potential 
energy and the respective gluelump mass. This relation was first quantified in Refs. 2|, |3| and 
recently exposed in an unified and model-independent framework {4, 5| usingthe effective 
field theory "potential nonrelativistic QCD" (pNRQCD) j^l (for a review see p]). 

The physics of the static hybrid system in the short distance limit is governed by, at 
least, two physical scales. One is the soft scale ~ 1/r, the inverse distance between quark and 
antiquark, the other is the ultrasoft scale ~ = Vq—Vs, where Vo and Vg are the static octet 
and singlet potential, respectively. The effective theory pNRQCD is particularly suitable to 
study this limit, as it profits from the large scale separation: 1/r ^ AV. For small distance 
r the static hybrid state can be understood to consist of a color octet quark- ant iquark state 
(acting as a static source in the octet representation) coupled to some ultrasoft and/or non- 
perturbative gluons to form an overall color singlet state. A close relative of the static hybrid 
(as we will see below) is the gluelump, which consists of a static octet source (of whatever 
origin) attached to nonperturbative gluons forming a color-singlet state. 

In order to gain a deeper understanding of the physics of the static hybrid system one 
can also consider how it is qualitatively affected by changing the number of dimensions from 
four (4D) to three (3D) or two (2D). The three-dimensional result is also important on 
its own. Four dimensional thermal QCD effectively undergoes a dimensional reduction for 
large temperatures. Therefore, determining the renormalization group (RG) structure for 
the static potential in three dimensions might open the way to a resummation of logarithms 
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at finite temperature. Three dimensional space-time is moreover a good testing ground for 
renormalon issues, since the linear power divergences associated with renormalons in four 
dimensions become logarithmic divergences in three dimensions and can be traced back using 
dimensional regularization (see Ref. Last but not least the computations in less than 
four dimensions represent important consistency checks of the theoretical framework used to 
describe the 4D hybrid (and singlet) systems at short distances. 

The color-octet static potential is obtained by integrating out the soft scale and has re- 
ceived quite some attention in the last decade. For D = A, it was computed with two-loop 
precision in Ref. 8], the leading three-loop logarithms were obtained in Ref. ol, and the 
associated resumation of logarithms was carried out in Ref. lOj. The 4D two-loop result 



of Ref. Ml was confirmed in Ref. 



Ill ] by computing the correlator of two Polyakov loops at 



finite temperature, which is then related to the singlet and octet static potentials. From the 
terms given in Ref. js] we have also been able to read off the general D dimensional expres- 



sion for the static octet potential, and verified it using the method developed in Ref. 
This expression is infrared divergent for D = 3. The divergences must be canceled by the 
ultraviolet divergences of the ultrasoft contribution. We calculate the latter in pNRQCD for 
arbitrary D up to two loops and present the results in full detail on a diagram-by-diagram 
basis. We also compute the D dimensional ultrasoft two-loop correction to the static singlet 
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13|. The 



energy providing an independent confirmation of the results obtained in Refs. 
one-loop ultrasoft singlet and octet contributions are basically identical except for trivial 
modifications. The ultrasoft two-loop octet computation has previously been considered in 
Ref. 1^, where a one-to-one correspondence to the singlet result in the temporal {Aq = 0) 
gauge was claimed. This statement can be rather problematic due to the singular nature of 
the gluon propagator in this gauge. Actually, for D = 4, the anomalous dimension of the 
octet static potential we find differs from the one obtained using the relation between singlet 
and octet computation based on the Aq = gauge argument of Ref. [l^. In three dimen- 
sions the problem is even more acute, as the result obtained from the suggested singlet-octet 
correspondence does not have the correct divergence structure to make the computation 
renormalizable and to derive a meaningful anomalous dimension. Our explicit ultrasoft two- 
loop result in Feynman gauge, in contrast, does, and exactly cancels the infrared divergences 
of the static octet potential for D = 3. This represents a strong cross-check of our calculation. 
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The discussion in the present paper closely follows the lines of Ref. [ij. Therefore, we will 
skip some details by referring to that reference. The outline of the paper is as follows. In 
Sec. ini we introduce the theoretical setup and present the bare two-loop soft and ultrasoft 
results for the color octet. The expressions for the relevant pNRQCD diagrams in Feynman 
gauge are shown in the Appendix, both for the color octet and singlet case. In Sec. IIIII we 
obtain the static hybrid energy at next-to-next-to-next-to-leading logarithmic (N^LL) order 
(apart from the soft three- loop matching constant) for D = 4. Sec. IIVI contains the complete 
RG improved ultrasoft calculation up to NNLL level in three dimensions. In Sec. |V] we 
discuss the situation in two dimensions, and in Sec. IVH we present our conclusions. 



II. PNRQCD 

Up to the next-to-leading order (NLO) in the multipole expansion (and irrespectively of 
the space-time dimension) the effective Lagrangian density of pNRQCD in the static limit 
takes the form^: 

C = Trjst {tdo - Vs{r)) S + {iDo - K(r)) O 

+ ^^A(r)Tr{OV-ES + SV-EO} +^^5^^Tr{0Ur-E,0}} + 0(r2). (2) 

We define color singlet and octet fields for the quark- ant iquark system hj S = S{r, R, t) and 
Qa _ Qci^Y^ j^^ respectively. R = (xi + X2)/2 is the center position of the system. In order 
for 5* and to have the proper free-field normalization in color space they are related to 
the fields in Eq. ([2]) as follows: 

1 

S = ^^, = ^=0", (3) 

where Tp = 1/2 for the fundamental representation of SU{Nc). All gluon fields in Eq. ([2]) 
are evaluated in R and the time t, in particular the chromoelectric field E = E(R, t) and 
the ultrasoft covariant derivative z-DqO = idoO — g[Ao{Il,t), O]. 



^ In principle, to fully account for the nonperturbative contributions at ©(r^) we have to add the term 



S£ = gY£^r\^Ti {O^ [D*E^ O] } (1) 



to the Lagrangian in Eq. ([2]). 



In the following we will use the index "S" to explicitly denote bare quantities. Parameters 
without this index are understood to be renormalized. We will furthermore use the notation 
D = 4: + 2e = n + 2e„, where e„ = -^^j^ parameterizes the (typically infinitesimal) difference 
to the closest integer dimension n = 4, 3, 2. The bare parameters of the theory are the 
coupling constant ob (qb) and the potentials V{s,o,A,B},B{r), generically denoted by Vb- The 
associated renormalized coefficients and V{s,o,A.B}if'] i^) are the Wilson coefficients of 
the effective Lagrangian and depend on the renormalization scale (z/). They are typically 
fixed at a scale smaller than (or similar to) 1/r and larger than the ultrasoft and any other 
scale in the problem by matching the effective to the underlying theory, which in this case 
is QCD in the static limit. 

In our convention has integer mass dimension, [a^] = [a] = M^~" {[9b] = M^"-^), 
and is related to (7^ by 

«B = ^, (4) 

where u is the renormalization scale. It has a special status since it does not receive correc- 
tions from other Wilson coefficients of the effective theory. Therefore, it can be renormalized 
multiplicatively: 

Ob = Zaa , (5) 



where 



The RG equation of a is 



In the limit — > 



s=l 



z/— a = a/3{a; e^) = 2e„Q! + a(3{a; 0) . (7) 



v^a = a/3{a; 0) = a/3{a) = -2a-^zW . (8) 
du da 



For the octet potential we employ an additive renormalization convention: 

Vo,B = Vo + 6V0 . (9) 

The counterterm 6V0 generally depends on the Wilson coefficients of the effective theory, i.e. 
on a and V, and on the number of space-time dimensions. Using the minimal subtraction 
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(MS) renormalization scheme in D{n) dimensions we define 

oo ^ 

^^o = E4?^- (10) 



s=l 



From the scale independence of the octet bare potential 



iy-^Vo,B = 0, (11) 



one obtains the RG equation of Vq. It can schematically be written as: 

u^Vo = B{V), (12) 

B{V) ^ - (i^-^SVo^ ■ (13) 

Note that B(y) is, in general, a function of all the potentials appearing in the Lagrangian. 
Note as well that Eq. ( fT2l) implies that all the 1 /e„ poles disappear once the derivative with 
respect to the renormalization scale is performed. This imposes some constraints on 6Vo: 



0(l/e„) : B{V) = -2«^Z« , (14) 

0{l/el) : B{V)^Z^] + «/3(a)|-Z« + 2«|-4f = , (15) 



and so on. 



A. Bare results in D dimensions 

This subsection summarizes the bare results of the calculations relevant in this work. In 
D dimensional momentum space we can write the bare octet potential as 

^ n=0 ^ ' 

where c"q\d) = 1 and c^i\D) = c^{'\d), c^i\D) being the analogously defined one-loop 
coefficient for the static singlet potential, see Ref. [isl. The two-loop coefficient c^\d) 
differs from c^\d), which was computed in Ref. [l^. We denote this difference by 

6c2{D) = i''\D) - di'\D) . (17) 
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ll| using thermal 



^02(4) was first obtained in Ref. jsj. The result was confirmed in Ref. 
QCD and effective theory methods to compute the correlator of two Polyakov loops and 
relate it to the singlet and octet potentials. Although no explicit expression is given one can 
also read the D dimensional coefficient Sc2{D) off the results in Ref. jsl. Using the approach 



developed in Ref. UJ] we have been able to confirm that expression for arbitrary D. It reads 



6c^\D) = Cl- 



7r-^r(e+i) sec2(7re) (e(2e + 3)r(e-i)r(2e+l) csc(7re) + 7rr(3e + |) sec(7re)) 



r(|-e)r(2e+l)2r (3e + |) 



After Fourier transformation to position space Eq. f lT6|) becomes (see e.g. Ref. 17l |) 

1 



2N, 
1 



2 ^ ^|v-2(n+i)e g(o)p) Y[l/2 + {n + l)e 
-9b — 

n=0 



■9l 



r (47r)"«/2 22-2ne^3/2+er[l _ ne] 

n=0 



(1^ 



(19) 



We now consider the ultrasoft bare correction to the hybrid energy at (9(r^) in the mul- 
tipole expansion. It can be determined from the calculation of the bare octet self-energy 
Ti^{E) in terms of (IPI) loop diagrams in pNRQCD. The octet propagator including the 
ultrasoft corrections (but neglecting nonperturbative ones) then takes the form 



dte'''\0\t)0\Q)) 



(20) 



To extract the ultrasoft correction to the static octet energy^ at (9(r^) it is sufficient to set 
E = Vg^\o{rO) in the bare expression for the self-energy^, i.e. SE'^% = ^ YT^jE = V^) + 0{r^^ 



The one-loop result equals the analogous singlet correction 
of the color factor and the replacement = Vo — Vg ^ — 
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19| except for a change 



loop) 



-9 



2Nr. 



-V(l 



r[2 + e]r[-3-2e] 



vr 



2+e 



(21) 



^ This correction is understood to be complex as long as we do not explicitly disentangle the real part associ- 
ated with the physical energy and the imaginary part associated with the decay width, see subsection III B I 
In order to consistently treat the 3D infrared divergence at ©(r") (which is canceled by an ultraviolet 
divergence of nonperturbative origin) in dimensional regularization it is important to set E = before 
performing the loop integrations of the self-energy diagrams in this "bare" approach. As a check we 
have computed the corresponding ofF-shell diagrams and explicitly performed the renormalization of the 
potentials and the octet field. Finally taking the "renormalized" on-shell limit E ^ Vo gives the same 
result for the ultrasoft correction to the static hybrid energy. 
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^ ^ 



FIG. 1: One-loop octet self-energy diagram at 0{r'^). 

In Feynman gauge it comes from the pNRQCD diagram shown in Fig. [H The two-loop 
expression is new and reads 

6E-{2 - loop) = g'^C^V^n-^ - 4e]^^r^[l + e] (22) 



X [^(e) - (1 + 2e)gr{e) + ^,(e)] {-AVb 



.2 AT/_^3+4e 



where 



^ 26^ + 66^ + 86 + 3 _ 26r(-26-2)r(-26-l) 
^^^^ 6 (262 + 56 + 3) (26 + 3)r(-46-3) ' ^ ^ 

^ 66^ + 1762 + 186 + 6 4(6 + l)nyr; 2 (6^ + 6 + 1) r(-26 - 2)r(-26 - 1) 
^'^^^ 62(262 + 56 + 3) 6(26 + 3)iV, 6(26 + 3)r(-46 - 3) ' ^ ' 

2-^-^(26 + 3)(262 + 26 + i)r(-26 - |)r2(-6 - |)r(26 + 1) 
^"(^^ = 6r(-46-3)r2(6 + i) • ^25) 

It is the central result of this work. In the Appendix we separately present the analytic 
expressions for the two-loop self-energy diagrams involved in the calculation of Eq. (!22|) in 
the Feynman gauge. We have also computed the corresponding two-loop expressions for the 



color singlet, which provide an explicit confirmation of the result obtained in Refs. 12|, [13 1 
and can also be found in the Appendix. 

The computation of the ultrasoft octet self-energy has also been addressed in Ref. jl^ . 
where it was argued that the result could be related to the singlet one (which is gauge invari- 
ant) in the temporal {Aq = 0) gauge. In order to reach this conclusion, the effects involving 
the gluonic content of the hybrid in the soft and ultrasoft computation were neglected (see 
Eq. 17 in Ref. [l4|). Then, after choosing the Aq = gauge, it is indeed easy to see that 
all non-vanishing self-energy diagrams at 0{r^) are, apart from the color factor and the sign 
of AV, the same for singlet and octet. Requiring gauge invariance of the final outcome one 
thus would obtain 

rp /•OO 

^^ {D-l){m-l)N/ I {vac\gEUmt, 0)^^ gE'^mvac) (26) 
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for the bare ultrasoft contribution to the static octet energy (in Euchdean spacetime) at 
(9(r^) and to all orders in a. 

At 0{a) (one-loop) Eq. ([26]) agrees with Eq. (EH). However, at 0{o?) (two-loop) Eqs. ([22]) 
and ( l26l) are different, since 5'o(e) = in Eq. ( 126]) . This is in contradiction with our non-zero 
result for 5'o(e) in Eq- (125]) . which was obtained by an explicit two-loop calculation in Feynman 
gauge (see App. lA 1|) . The existence of a non-zero value for go is required on a theoretical 
basis in Sec. IIVI Otherwise SE^"^^ would not produce the right divergence structure in D = 3 
dimensions to make the theory renormalizable nor lead to a consistent RG equation for the 
potential, as there would be some divergences proportional to logarithms of the ultrasoft 
scale left after the subtraction of the one-loop subdivergences. Our explicit computation, 
in contrast, is consistent with the renormalizability of the theory, and moreover perfectly 
cancels the infrared divergences of the soft computation. This represents a strong check of 
our computation. In Sec. IIIII we will see, that with our result in Eq. ( 122]) we also obtain a 
different NLL anomalous dimension of the octet potential Vo in four dimensions compared to 
Ref . [l^ . The corresponding cancelation of soft infrared and ultrasoft ultraviolet divergences 
in four dimensions can however not be checked, since the required soft four-loop result is not 
known at present. 

Actually, the fact that there are problems with the naive usage of Aq = gauge becomes 
apparent at several points of the calculation. In pNRQCD this already happens at 0{r^) 
in the multipole expansion for D = 3. At this order the bare soft result for Vo is infrared 
divergent (cf. Eq. (147])). but in the Aq = gauge there are no ultrasoft octet self-energy 
diagrams to cancel that divergence. It is also remarkable that in this gauge the space-like 
strings of the rectangular Wilson loop contribute to the static potential, whereas in Feynman 
gauge they are expected to be negligible^. This can be easily visualized for D = 2, since then 
time and space coordinates can be symmetrically interchanged (in the Euclidean), showing 
that in the Aq = gauge the contribution to the singlet potential comes exclusively from the 



space-like strings. Actual 



ythat is true for arbitrary D at tree level. For a detailed analysis 



at one loop see e.g. Ref. 20[ and the references therein. Overall, we think it would be very 



interesting to quantitatively study the effect of the asymptotic gluonic degrees of freedom 



In Ref. 



15| it was explicitly shown that the space-like strings in Feynman gauge do not contribute to the 



4D static singlet potential through two loops. 
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on the static singlet and hybrid energy in the Aq = gauge. 

To the above discussion one should also add the (9(r^) ultrasoft contributions from octet 
self-energy diagrams without singlet-octet transitions. They are proportional to or Vc 
and may have logarithmic ultraviolet divergences, relevant for the determination of the coun- 
terterm {6Vo) and the anomalous dimension of Vq. These divergences can, however, not be 
proportional to positive powers of AV^, because the ultrasoft scale AV^ does not appear in 
the diagrams with onshell external fields, and the corresponding contributions vanish in di- 
mensional regularization. Therefore they cannot solve the Aq = gauge problem described 
above, since the first discrepancies between our two-loop result in Eq. ( l22l) and the 0{a^) 
term in Eq. ( I26l) are linear in A^ for D = 3 and cubic for D = 4. 



B. Observables 



At short distances we approximate the bare propagator of a generic hybrid system as 
follows: 



dte-'(H'O"(t)H^O\0)} ~ ^_,,,_'^^(^^^,^ . (27) 

if" stands for the gluonic content of the hybrid {H°- = B"',E"', etc.) and T,b{E) is the 
self-energy of the system expressed in terms of bare parameters (potentials, couplings). It 
accounts for effects at the ultrasoft and the nonperturbative scale. 

Let us now construct some observables. In this paper we focus on the hybrid static 
energy and decay width. Their definition can become gauge dependent for unstable particles. 
Fortunately this is not so in our case, since both the on-shell and pole definition (see e.g. 



Ref. (2]J for the respective definitions) are identical at the precision of our computation. 
Close to the resonance we can expand T,b{E) around the complex pole position Eh + iTH/'2 
of the hybrid propagator in Eq. flTTI) . In this paper we will consider ultrasoft self-energy 
contributions S^* up to two-loop level and NLO in the multipole expansion. It is therefore 
sufficient to approximate 12b{Eh + ^r_f//2) by Sb(£'°) and set Eq = V^\ci{a2) for D = 4 and 
Eq = V^\o{a) for D = 3. Note that in the latter case C(a^) terms would produce subleading 
0{r^) corrections to the energy and decay width. 

We now choose to shift the energy origin by adding —2mB, i.e. E ^ E — 2mB- The 
term 2mB = 2m -|- 26m(^soft) will allow to absorb ultraviolet soft divergences (into 6m(^soft)) 
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that may appear in the computation, m has a clear physical meaning in the case of heavy 
quarkonium: it represents the heavy quark (pole) mass. At NLO in the multipole expansion 
we thus have 

g _ - V^B \ - 2^,) - E - E^l .r„/2 ^0((E-E„^ ,:r„/2)°) . (28) 

where is the normalization of the hybrid propagator and 

EH-iTH/2 = mB + V,'' + Eb{Eo) (29) 

is the complex position of the hybrid propagator pole. It can be understood to be composed 
of contributions from different energy regions. After factorization we can express it in the 
following way (either bare or renormalized) : 



Enir) - iTH/2 = 2mB + Vo,b + SE^% + 5Ej 



B 



The subscript MS of ?71ms(^s) denotes the scheme for possible hard infrared divergence sub- 
tractions (canceling the ultraviolet soft divergences), e.g. in three dimensions, and should 
not be confused with the ultraviolet (hard) renormalization scheme for the mass. In four 
dimensions we have rriyi^^h's) = mos, where rrios is the usual on-shell mass of the heavy 
quarks. Similarly V^jj^^r; z/^, z/) — )■ V^jj^{r] u) in four dimensions, where is the octet po- 
tential and encodes the effects due to soft degrees of freedom. SE"^^ and 5E^ encode the 
effects due to the ultrasoft and nonperturbative degrees of freedom respectively. Each term 
depends on the respective factorization/renormalization scales: Vs separates the hard (m) 
from the soft scale (1/r), v separates the soft from the ultrasoft scale (A\^) and z/„p separates 
the ultrasoft from the nonperturbative scale (Aqcd)- For real factorization scales Vs-, i^np, 
the decay width F^^ is contained in the terms 6E^^ and 6E^. We give more details in the 
next sections. 

III. RESULTS FOR D = A 

In four dimensions the renormalized coupling constant a satisfies a non-trivial RG equa- 
tion. We will work here in the MS renormalization scheme, which is related to the MS scheme 
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by a redefinition of tlie renormalization scale: z^(ms) ~^ ^(ms) "^m^' wliere c 



-MS 



62 



l(ln(47r)-7B) 



Tlie counterterm coefficient 



is understood in the MS sclieme, but tlie form of tlie anomalous dimension of a, 

„2 



(31) 



Q (X 

a/3(a) = -2a(/3o— + 



+ . . . 



(32) 



is the same in both schemes, MS and MS. The constants /3o = f^F^-f and f3i 



34 r^2 



C\ — ACpTpUf — ^CXTrnf are the standard MS (MS) one- and two- loop coefficients of 
the QCD beta function. 

We also have to consider possible soft and ultrasoft corrections to Va- They were obtained 

n n n 

in Ref. [10[ with LL accuracy, in Ref. (13| with NLO accuracy and in Ref. [iJ] with NLL 



accuracy. The outcome is 



p^Va = + 0{a^) 
dv 



(33) 



for the anomalous dimension and Va = 1 + 0{a^) for the initial matching condition. We 
conclude that for the precision of our calculation we can use = 1. 

The counterterms of the octet potential, which subtract the ultraviolet divergences from 
the result of our ultrasoft two-loop computation in Eq. fl22l) in four dimensions, read 



7(1) 



7(2) 



a_ [Ca(47 - 127r^) - WTpUf] 
3^ ^ 



IOStt^ 



2 
3 



a 



(34) 
(35) 



2A^/"' 3"^^ (47r)2 ■ 

The latter expression comes from the l/e| pole of the two-loop result and from the a'^/e^ 
divergence of as in the divergent term of the one-loop self-energy. 

From Eqs. ( IMl) and (l35l) we can derive the RG equation of the octet potential at two-loop 
order^. We find 



(36) 



In our counting this translates to N'^LL order, because AV comes at least with one power of a. 
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This result holds in any momentum-independent renormalization scheme that is related to 
the MS scheme by a {D independent) redefinition of the renormalization scale u. Solving 
the RG equation we can write the RG improved static octet potential as 



where 



X <^ — In— — - (a(z/) - Q!(i^i)) ' 



37r a{Pi) ' ' V3(4vr)2/3o 277r2 J 

describes the ultrasoft RG evolution of and the initial matching condition at the (soft) 
scale z/j is given by 

v:,feen.)-5^^i:(tr)"°^°'(^^-) 

n=0 ^ ^ 

with coefficients (aQ°''(r; z/j) = 1) 

a(\r; Ui) = ai + 2/3o In {uie""^r) , 

2 

4°) (r ; u,) = ai"^ + ^(3,' + ( Aa,(3o + 2/3,) In (z/,e^-r) + A/3^ In' (u.e^^r) , 

+ (^27i'/3^ + 6ai"^/3o + 4a,/3, + 2/^2 + y C'jTr^^ In (z/^e^^r) 

+ (^12ai/3o' + 10/3o/3i^ In' (z/,e^^r) + S/3^ In^ {v.e^^r) . (40) 



Explicit expressions for a, and can be found in the literature 16|, |22|. The constant 



a~^^ is still unknown and represents the only missing piece in the N^LL result for the Wilson 
coefficient Vo, Eq. ( l37l) . Note that this expression is real for positive renormalization and 
matching scales v and z/j. 

We can now determine the complex pole of the hybrid (octet) propagator in four dimen- 
sions. Up to nonperturbative effects, which we neglect here, it is given by 

r r 

^uir) - ~ Eo{r) - - 



13 



Therefore, besides the octet potential, we also need the ultrasoft correction in the MS scheme, 
which we obtain from our bare two-loop result in Eq. (122|) after subtraction of the divergences. 
It reads 

. . 1 



2Nr 





'-AV 


+ 61n2-5^ 


"-^fein 


97r V 


V 





(42) 



+ 



a- 



1087r2 







(l8/3oln' 


V 



-6(iVc (IS-Stt^) -2/3o(-5 + 31n2))ln 



V 



- 2iV, (-84 + 39 In 2 - 87r2(-2 + 3 In 2) + 72C(3)) + /3o (67 + 'iix'^ - 60 In 2 + 18 In^ 2) 

Note that this object is complex for positive v. Since the octet potential is real, the imaginary 
part of (5-E'"|jg(z/) directly gives the decay width. Choosing v = AV in Eq. fj42l) we find for 
the decay width of the hybrid (octet) system 



3 ^ '2Nr 



1 + ^it^ {Nc (13 - 871^) - 2/3o(31n2 - 5)) + 0{a') 
iZtx 



where we have to use the NLO expression for the ultrasoft scale AV: 

Ca 



AV{r- v>) 



n=0 ^ ' 



(43) 
(44) 



which is scheme and factorization scale independent at this order. 
For Eh we need the real part of SE^^^^u): 

Enir) = V^-Ur; 1/r) + SV%{r- 1/r, AV) + Re SE^^^^AV) . 



(45) 



Replacing the first two terms in this equation by Eq.( !371) and using the leading term of 
Eq. ( l42l) in the last term, we reach N^LL accuracy. Eq. ( l42i) also provides the subleading 
ultrasoft correction to the hybrid energy relevant at N^LL order. 

We will not attempt to compare our results with 4D lattice simulations. This would 
require the incorporation of nonperturbative effects and the proper treatment of renormalons, 
which goes beyond the scope of this paper. 



IV. RESULTS FOR D = 3 



The derivation of the static octet potential in three dimension is quite analogous to the 
color singlet case, which has been discussed in detail in Ref. We therefore focus on the 
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novel aspects of the analysis for the octet (hybrid) system in this section and refer to Ref. [ij] 
for a more careful account on the universal issues related to the 3D static potential. 
As argued in Ref. [ij], the coefficients Vaib are not renormalized at (9(r°), i.e. 

Zaib = 1 + 0{r) . (46) 

The reason is that the potentials and a have to appear perturbatively (with positive powers) 
in the counterterms, otherwise the renormalizability of the theory at leading order of the 
multipole expansion would be spoiled. Moreover, we can set Va = Vb = I (just like in 
4D), as 0{a) soft corrections would be multiplied by factors of r and would move us away 
from the precision of Vo aimed for at this paper. Actually, from inspection of the possible 
diagrams that will contribute at the soft scale, we know that Va/b = 1 + C^('^^) |l3 |. 

From the bare soft computation we can completely fix 6Vo through 0{r^). We obtain (for 
the ultrasoft counterterms) 

4? = ^ « + r^AV^a {C, - ^) i + r^AV CAa^ (c,-^)\ (47) 



-rV(C 



X 



_Ca\ 
2 ) 

137r2 - 2304) + 8 (lOvr^ + 144) C^T^rij - ASTpUj (4 (tt^ - 10) Cp + ir^TpUf) 

2304 ' 

Zg) = r^AVa^C, - ^jc^ + r^a^[c, - , (48) 



zg)=rV (49) 



TJ^^48 

= V n>3. (50) 

The soft calculation is organized in powers of ar. The tree level computation gives the 
first term in Eq. fHTl) (once the ultraviolet divergences of the soft one-loop heavy quark 
self-energies have been subtracted). The soft one- loop contribution to the octet potential is 
infrared safe. The two-loop computation leads to the remaining terms. These results are 
exact at C(r^). 

The fact that one can renormalize the potential with a finite number of terms at a given 
order in the multipole expansion, i.e. ZyJ = for n > 3 at C(r^), refiects the super- 
renormalizability of the theory in three dimensions. 
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FIG. 2: One-loop contribution to the octet propagator at 0{r^). The dotted line represents the 
field. 

The infrared divergences of the bare potential have to cancel the ultraviolet ones of the 
ultrasoft computation. We can then use the explicit one- and two-loop ultrasoft contributions 
to partially check the structure of the octet potential counterterm in Eqs. P71) - P9l) . 

At leading order in the multipole expansion the octet field has a residual interaction with 
the ultrasoft gluon field. The octet potential receives an ultraviolet divergent correction from 
the one-loop self-energy diagram shown in Fig. [21 which yields the first term in Eq. (147]) . 
This is the only term at 0{r^), since higher loop diagrams cannot contribute because the 
coupling and the potentials have to appear perturbatively (with positive powers) in the Z's. 
Since a has dimensions of mass, the potentials would appear with negative powers in multi- 
loop diagrams at 0{r^). This is not allowed by renormalizability. By the very same reason 
the octet field does not require renormalization at 0{r^): 



We now consider higher orders in the multipole expansion. At two soft loops in dimen- 
sional regularization infrared poles up to 0{l/el) appear. The ultrasoft computation in 
pNRQCD yields the following results for the counterterms: 

1) The second term in Eq. (147]) comes from the 1/e-s divergence of the ultrasoft one-loop 
correction in Eq. fl2Tl) . It is scheme independent and fixes, together with Eqs. f|T^ and f|T5|l . 
the first term of Eq. (148!) and Eq. (l49l) . It would also be possible to compute these l/e| 
and l/e| divergences directly. The l/e| term is confirmed by the ultrasoft two loop result 
in Eq. fl22p . The l/eg term would however require an ultrasoft three- loop calculation, which 
has not been performed yet. 

2) The third term in Eq. follows from the remaining l/e^ divergence in the ultra- 
soft two-loop correction, Eq. (l22l) . once all sub divergences (associated with the 0{r^) octet 
potential) have been subtracted. This result combined with Eq. flTSl) then fixes the second 
term of Eq. fjiSll. 



(51) 
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In summary: the explicit ultrasoft computation allows to check the three first terms of 
Eq. f H7|) and the first term of Eq. fHHl) . We find perfect agreement. The use of Eqs. f|T4|) 
and (fTSl) allows us to completely check Eq. ( l48l) and ( l49l) . We also find perfect agreement. 
Note that this can be understood as a non-trivial cross-check of two independent determi- 
nations of these terms. 

The only remaining term is the last one in Eq. (H7|) . This term would be canceled by 
the sum of the ultrasoft ultraviolet divergences from three-loop octet self-energy diagrams 
with and without singlet-octet transition. The latter diagrams are proportional to V| or Vq 
and scaleless (for onshell external fields), as they are insensitive to the ultrasoft scale AV. 
Therefore, the associated RG evolution will run down to the non-perturbative scale a. This 
is not so for the ultrasoft diagrams with singlet-octet transition, which get infrared regulated 
by the ultrasoft scale correspondingly producing logarithms of AV. Without an explicit 
ultrasoft three-loop calculation we are not able to distinguish among the divergences from 
the two types of diagrams nor check the last term of Eq. f l47|) . which has been deduced 
from the soft result. Therefore, we turn the problem around and use the latter to fix the 
ultraviolet divergences of the ultrasoft three-loop computation for the next sections. 



A. pNRQCD RG 

We can now deduce the RG equations of Vq. Again the discussion is similar to the one of 
Ref. to which we refer for extra details. 

At leading order in the multipole expansion the RG equation of the octet potential reads 

ly-^Vo = -Cacx . (52) 

The running of AV at this order reads (including the tree level matching condition) 

AVxir; v) = -uCa Hrvdx) + 0{e) , (53) 

where X stands for the factorization scheme (e.g. MS or MS): 

dms = e^^/^y^ ~ 2.36546, d^ = c/ms = e^^/2 ~ 0.890536, (54) 
withcTTo =ei/2(W^'^)-7B). 
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Now we consider the sub leading contributions to Vo- The complete anomalous dimension 
of the static octet potential through O(r^) takes the form 

2 

V: = R(V^ = -Cah + ' 

du 



u^Vo = B{Vo) = -CAa + CpaVyaB^{CAaT{AVf-^, (55) 

n=0 



and is explicitly given by 



v^V^ = -CAa - r'AVy[Cj. - ^)vf- - 2r'AVx CAa'(Cj. - ^)vl (56) 



+ r'aHC,-^W 



X 



(137r' - 2304) + 8 (IQtt^ + 144) CaTfUj - iSTpUf (4 (tt^ - 10) Cp + n^TpUf) 

384 



The form of the anomalous dimension in Eq. flS^ is invariant under scheme transforma- 
tions that amount to a redefinition of the renormalization scale: u — )■ i^c]^^, where cx is 
a ea-independent constant (cms = 1)- We indicate renormalized quantities in this class of 
momentum independent renormalization schemes by an index "X" (X=MS, MS, etc.) as 
e.g. AVx in the above equations. 

By solving the RG equations in the MS scheme we obtain the static octet potential with 
O^a^r"^) precision: 



Vo,Ms{r] Us, i^) = V^Msir; i^s, i^i) + ^V^Msij^', ^) , (57) 



where 



SVJSsir; r, ^) = - CAa ln(z/r) + {Cp - ^]rW 

1„.. . 1 



X < - -Cl ln\ru) - -Cl [ - 4 + + In vr] ln\ru) 



+ 



+ ^( - 7l + 87E - 48 + (8 - 27^)ln7r - InV) 



384 8 

+ n,T,(c^(3 + ^) + ^-(5 - y)) - (-/T.)^y 



ln(rz/) } (5^ 
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is the running and 
1. 



TT 25 31 

^(397E - 715 + 564 In 2 + SQlnvr) + - - -((3) 



- — (7^^ + In vr) (144 - 127^^ + 7I + In^ vr - 12 In vr + 27^^ In tt) 



TT^ 1 

— (15-67£-81n2-61n7r) + -(57B-8 + 51n7r) 



vr^ 1 
-n/TpC^ ;^(577ij - 97+ 12 In 2 + 571n7r) + -(97^ - 22 + 91n7r + 10C(3)) 

-2oo D 



+ ('^/^F)'[^(7-37i.-ln(167r3)) 



(59) 



is the initial matching condition. Note that the tree level and one-loop matching conditions 
have been included. The tree-level result depends on the factorization scale i^s ^ 1/r, which 
separates the hard and the soft regime. The dependence on this factorization scale cancels 
the infrared scale dependence of the mass, cf. Eq. fl30|) . 

We stress that we have obtained the exact contributions to the static potential through 
O(r^). There is nothing left. Moreover, by setting v ~ ISV large logarithms are resumed up 
to the scale A^, i.e. we have determined all logarithms of rAV" at O(a'^r^). 

The above results have been presented in the MS scheme. We can easily transform 
them to another momentum independent scheme by redefining the renormalization scale. In 
particular, if we write the expression in terms of Al^, most of the scheme dependence gets 
encapsulated in A^. For instance, Eq. (155]) can be reexpressed as 



1/CV _ 1 
6 VC^ 2 



r a 



Vl[^V^{r-v)-AV^{r-l/r)) 



+ -'^A«'(§; - (AVi(r; u) - AV^r; 1/r)) 



2 2f^F 1^ta2 



X 



;i37r2 - 2304) + 8 (197r2 + 144) C^Tpn/ - ASTpUf (4 (tt^ - 10) + vr^T^ri;) 

384 



X {AVx{r;u)-AVx{r;l/r)) 



(60) 
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B. Ultrasoft contributions {SE"^^) 



The general structure of the ultrasoft contribution SE'^'^ at (9(r^) is 

n n+1 



AV 



(61) 



n=0 ^ ^ s=0 

where we have set z/„,p = i^. The dependence on Unp first appears at three loops in this (9(r^) 
contribution to SE"^^. At present, concrete results for the ultrasoft corrections are available 
at one, 0{g^), and two loop level, 0{g^), given in Eqs. fl2T|) and fl22l) . After MS subtraction 
they read 



5E-MS(1-100P) 



a Vir' AV^s U + 7i? - ln(47r) + 2 In 



MS 



(62) 



5E:Ls(2-1oop) 



■a^V^r^ AVms 



CaIii' 



-AK 



MS 



+ CA(2-7i? + ln(47r))ln 

11 \ 

- 2CA(vr' -7+-{2-7E + ln(47r))2) - 2n/T^ j . 



(63) 



Note that the these expressions have imaginary parts. Setting u ~ AVx{i^) resums large 
logarithms into the potential Vo and minimizes them in the ultrasoft contribution SE'^'^. It 
also simplifies the determination of the decay width of the octet (hybrid) system: 

^ 1 ■.-^0 0»-.-^0/_\ 1 



H 



a V'i r'AV^^{v^s) vr + -—a' V'i r'AVMs{^us) C^7r(2-7s+ln(47r)) + ^(aV^AyO^ 



27V, 



27V, 



where 



AVx{vus) = CAaW{l/{CAadxr)) 



(64) 



(65) 



and W{z) is the Lambert function (here: X = MS). 

The precision of Eq. f l64p is limited by the existence of some unknown single logarithms 
that appear in the three- loop ultrasoft self-energy (see next subsection). The imaginary 
parts of those logarithms produce terms like ST ~ a^r'^. 
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C. Subleading ultrasoft and nonperturbative effects 



Even without an explicit computation we can obtain some information on the ultrasoft 
three-loop terms. At three loops we start to have contributions from diagrams with no 
singlet-octet vertices. They are proportional to or Vc and vanish in dimensional regu- 
larization (for onshell external fields), as they are insensitive to the scale AV. Therefore, 
any ultraviolet divergence proportional to V| or Vc should be canceled by an infrared one, 
which signals a sensitivity to the nonperturbative scale a. We quantify this statement with 
the following equation 

u^SE^iu, v) = -B{V) - z/-^ J , (66) 

where the RG structure of the nonperturbative term is the following 

u^dEZ' = CAa + B r^a^l + 0(a/AV)) . (67) 
du 

Solving this equation we obtain 

SE'ifiv^,) = C^a (In (^) + ch) + BaV In (^) + 0{a'r') , (68) 

where ch is a nonperturbative constant that depends on the specific hybrid (gluelump) we 
consider. Actually, the 0{r^) term is nothing but the gluelump mass 

AhM = CAa [In (^'^^+ch). (69) 

It can be related to a gauge invariant correlator as 

Ah= lim Un{H'^{T/2)<p:l\T/2,-T/2)H\-T/2)) , (70) 

1 — ^oo ± 

where H"^ is the field operator associated with the gluonic component^of the hybrid and 

J,adj 



(T/2, —T/2) is a Wilson line in the adjoint representation, see Ref. 4|. 



Note on the other hand that B is independent of the hybrid type and can be obtained 
from perturbation theory. It is however unknown at present and constrains, besides the 
nonperturbative gluelump mass Ah, the precision of our result. 

At this point we could also analyse subleading ultrasoft effects in the a/AV expansion 
along the lines of Ref. as some of them, namely the logarithmic terms proportional to V|, 
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are fixed by the RG structure. In view of the dominant uncertainties discussed above, and 
because there are more unknowns than in the singlet case we refrain from performing 
that analysis here. 

Combining Eqs. ( 157|) . ( I62l) . ( 163|) and ( l68l) we can write down the static hybrid energy 
Eh = 2mMs+K,MS+<^-E"Mg+5-E'^^]y[g. Expressed as a double expansion in ar and 1/ ln(rA\^) 
it reads with O^a^r"^) and NNLL accuracy 

Eh = 2mMs{'^s) + Cpa ln{r^i^^^ne^^) + ^ (Cf - -f) (^Ca - infTp)a\ 
+ SV,%{r- ^, z7„,) + 6E:^usi^^s, i^np) + SE'if^MsM + 0{a\'' ln°) 



Ca 



MS 



\{CF-^){^CA-AnjT,)ah 

(Cf - ^)a\^\^Cl\n\rAVus) + \cI{2^e - 3 - 2 ln2) ln2(rA\4 

Cl f ^ + ^ (47I + 4 In^ 2 - 7^(10 + In 256) + 38 + 8 In 2 - 2 In vr) 
V 384 8 

- n;T^(c^(2 + ^) + C^(5 - y)) + (%T^) V ln(rA\/MSj 

+ 0(aVln(A\/Ms/«))- (71) 

We have checked that the explicit scheme dependence of ^V/^g and ^-E" ms ^^'^ ^'^^ implicit 
scheme dependence of Eq. f lTT]) through the logarithms of AV cancel up to (9(a^r^ \vP). Note 
also that cj^^ms is scheme dependent. 

The leading uncertainty of Eq. (17T]) comes from the nonperturbative gluelump mass Kh- 
The constant h^n is independent on the distance r and therefore drops out in the force, i.e. 
the derivative of the potential. In that case the leading uncertainty is due to the coefficient 
B, which is unknown but, unlike Kh-, of perturbative origin and could be determined by an 
ultrasoft three loop calculation. 

Eq. (17T]1 represents one of the main results of this paper. One could try to perform some 
comparison with the (quenched) lattice data existing in the literature 23|-|27|. It is not clear 
though that these simulations reach short enough distances so that our results can be tested 
quantitatively (on the other hand it should not be too costly to perform a dedicated short 
distance simulation to test our expression). Moreover, it would be as well very interesting 
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to study the accidental symmetries and degeneracies among different hybrid energies that 
should appear at short distances along the lines of the studies in four dimensions performed 
in Refs. jj, IsJ. To this end one could consider energy differences where the perturbative 
terms (the non-analytic, i.e. logarithmic dependence in r) cancels. The non-canceled terms 
should be produced by the nonperturbative terms (note that the scheme dependence of ch/w 
cancels in the difference): 

Eh - Eh' = CAa{cH - c^O + Cnwr' + 0{r') . (72) 
All this would require a dedicated analysis, which will be carried out elsewhere. 



V. RESULTS FOR D = 2 



It is also interesting to investigate static hybrid systems in two space-time dimensions. In 
the following we set = for simplicity. In exactly two dimensions physical hybrids do not 
exist, as there are no propagating (physical) gluons. It is instructive to see how this finding 
arises in an explicit calculation. For a typical hybrid, the object to be computed is e.g. 

(Wh) = ( ^Tr V E(T/2, R) ■ E(-T/2, R') exp {'wf^ A,dx^^) ) , (73) 

where the contour F of the integral in Eq. (1731) is a rectangle with spatial extension r and 
temporal extension T in Minkowski space and we have inserted two chromoelectric fields on 
each space-like end-string of the rectangular Wilson loop. We could have also chosen other 
gluonic configurations instead, but this would only make the discussion more complicated 
without changing the physical outcome. 

Eq. (!73|) is gauge invariant. Choosing axial {Ai = 0) gauge in exactly two dimensions, 
the only non- vanishing component of the gluon field-strength tensor is Fqi = —Fiq = —OiAq. 
Hence, [Aq) gluons neither interact among themselves, nor propagate in time, since no time 
derivative acting on the gluon field is left in the Lagrangian. Therefore only^ planar "ladder" 
diagrams (with "potential" gluons) contribute to (Wh)- In the large T limit the final result 

® We do not include in this discussion possible interactions of E with the gluons in the time-like strings. 
They would produce terms proportional to the singlet: ~ e*^"-^, but not contribute to the hybrid energy. 
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takes the form 



{Wh) ~ {E''{T/2)E\-T/2))e 



-iVoT 



(74) 



where V^(r) = i^ar. Note that this octet potential is the same as obtained from a leading 
order, tree- level, computation in 2D static QCD. Therefore the non-existence of the hybrid 
does not arise from the fact that the octet potential diverges in two dimensions (generating 
an infinite hybrid mass) but rather from the fact that 



is local in time. That is because the gluon does not propagate in time. 

For D > 2 one expects (Wh) to converge towards the previous result in Eqs. fTMl) . fl75|) 
as Z) — )■ 2, provided the limit is smooth. The discussion parallels to some extent the one 

n 

for the singlet potential in Ref. [![. Like in that case we cannot make definite statements, 
since the calculation of the Wilson loop becomes intrinsically nonperturbative. The reason 
is that the proper ultrasoft expansion parameter is now AV/g ~ gr, see Ref. Thus the 
nonperturbative ultrasoft (AV/g) expansion and the perturbative multipole (gr) expansion 
mix. 

All we can say is that the cancelation of the soft and ultrasoft one-loop contributions 

n 

we observed for the singlet potential [JJ also occurs for the hybrid potential, because the 
computation differs only in the prefactor. For the chromoelectric correlator we find at leading 
order in a 



This is what we would expect from the decoupling of gluons in two dimensions. Nevertheless 
at NLO in a one obtains a finite contribution. Again we can not draw any definite conclusion, 
as this result is obtained in standard perturbation theory, whereas the correct expansion 
entails inverse powers of the coupling constant g. 

VI. CONCLUSIONS 

We have computed the energy of a static hybrid in the weak coupling limit in four and 
three space-time dimensions. Employing finite temparature theory methods we have checked 




(75) 




(76) 
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the static potential of a color octet quark-antiquark pair with distance r at two loops for an 
arbitrary number of dimensions D. The result represents the soft contribution to the static 
hybrid energy. Using the effective theory pNRQCD we have explicitly calculated the missing 
ultrasoft two-loop contributions. We have also confirmed the respective two-loop result for 
the color singlet. The cancelation of soft infrared and ultrasoft ultraviolet divergences is a 
strong cross-check of our results. 

For £) = 4 we have determined the static octet potential and the static hybrid energy 
through N'^LL order, up to the unknown three-loop soft matching condition. The result 
disagrees with an earlier result obtained in Ref. We have also given the ultrasoft 

contribution to the static hybrid energy at N^LL and computed the decay width of the 
hybrid/octet system at NLL order. 

In £) = 3 dimensions our result for the static hybrid energy reaches 0{a'^r'^) accuracy 
in the (soft) multipole expansion. At this order we have determined the complete expres- 
sion at NNLL order in the ultrasoft aj ISV expansion. The precision of the result is only 
limited by unknown terms of order Oioc'r^Xwi^ISV j ay) and the nonperturbative gluelump 
mass. The former can in principle be obtained by a perturbative three-loop computation, 
whereas the latter requires lattice simulations (note, however, that the gluelump mass van- 
ishes in derivatives of the potential, like the force). Besides the energy we have determined 
the 3D hybrid/octet decay width through 0{o?r^^ in the multipole expansion and NLO in 
the aj ISV expansion. 

We have also studied the two dimensional case, where the exact result is known: hybrids 
do not exist. This is due to the temporal locality of the gluonic correlators as there are no 
propagating (physical) gluons in exactly two space-time dimensions, see Eqs. ( 174|) and ( 175|) . 
rather than due to an infinite octet potential or hybrid mass. In the D ^ 2 limit this result 
is more complicated to obtain, because already the ultrasoft contribution is intrinsically 
nonperturbative. Strong cancelations have to occur among the different contributions from 
the soft, the ultrasoft and the nonperturbative {g) scale. We can only report partial and 
inconclusive results on this limit, particularly because they are based on perturbation theory, 
which is not reliable for D — 2. Nevertheless, it might be worth exploring the D ^ 2 limit 
in more detail, as it could provide nontrivial information about the dependence of the {D 
dimensional) perturbative results on [D — 2). 
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Appendix A: Two-loop pNRQCD self-energy diagrams in Feynman gauge 



1. Color octet 




d\2 



r(5-2rf)r(f) 



° 4(rf-4)(rf-3) 



(Al) 



2L+1 



48 



1 + ^(ea) 



(4) 



1 I -6L+ll-61n2 _ 2L^ _ 2L(6 In2-ll) 
12iJ 18e4 3 9 



349 
108 



2 In^ 2 I 22 In 2 



+ 



9 




r(5-2d)r(f -1)' 

16 



(A2) 



a 



(3) 



_l__L+l + 0(63) 



Co 



(4) 



L I 11-6 ln2 

24e4 6 36 




(rf-2)r(5-2rf)r (f - 1)' 



1 

4^ 



8(ci-3) 

1-2L + 2L - ii?^ - 2 + C(e3) 



(A3) 



2e3 



24 



1 

6e4 



2L I 25 2 In 2 , /'/-)/', N 
~ + 18 ~ ~ + '-^(.^4j 




2i-2d (^2 _ ^ r (I - d) r(i - dfv [d + 1) v{df 
(rf-4)r(^)r(^) 



+ 



v^2 



-d-2 



-3rf^ + 30c/3 - 111^2 + 181^ _ 112) r(5 - 2d)T (f - 2) V{d - 4) 



(3) 



_L _ 2L + I (15 - vr^) + ^(ea) 



(A4) 



(4) 



_ + -9^-2.;+15-91n2 - L'' - \L (-30 + 471^ + 18 111 2) 

-|(8C(3) - 27 + 20 In 2 - 6 In^ 2) - ^£(192 In 2 - 47) + C(e4) 
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3(3 - d)T{5 - 2d)T (1 - f ) r (f - 1) r (f ) 

ler (3 - 1) 



a 



(3) 



3 + 0(e3) 



^ 3L-4+3In2 ^ ^2 + 2L(31n2-4) _^ 3^ ^ 26 ^ ^^2 ^ _ ^ 



6e4 




CA(3rf - 2) - 4(rf - 2)n/Tp 



X 



0F(d - 3)r(5 - 2d)T (1 - f ) T{d - 2)r (f ) 

2^+3r (3 - f ) r {^) 



a 



(3) 



I - ^ + ^(^3) 



In 2 
9 



. +[^-i](^^ + M21n2-|) + ^ + i+ln^2--|^)+0(e4 



L :=ln 



Ci'^ := - 



7r^ 
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2. Relation between singlet and octet diagrams 







AV^-AV 



CF^i^CA — Cp) 

AV^-AV 




Cf^(^Ca—Cf) 
AV^-AV 






(All) 
(A12) 

(A13) 
(A14) 



4-'^(rf2-6rf + io)r(i-d)r(i-f)^r(f + i)^r(rf+i) 



A. 




(rf-4)r(i^)r(^) 
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6e4 



+ ^ + V(31n2-2) + 0(e4) 



(A15) 



(A16) 





CF^i^CA — CF) 

AV^-AV 



CF^i-^CA — CF) 

AV^-AV 





(A17) 



(A18) 
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